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^ ■ Abstract: It has been shown that the string- flip potential model reproduces 

CN , most of the bulk properties of nuclear matter, with the exception of nuclear 

S-J \ binding. Furthermore, it was postulated that this model with the inclusion of 

^SJ ' the colour-hyperfine interaction should produce binding. In some recent work 

^^ . a modified version of the string-flip potential model was developed, called the 

flux-bubble model, which would allow for the addition of perturbative QCD 
interactions. In attempts to construct a simple qq nucleon system using the 
flux-bubble model (which only included colour-Coulomb interactions) difficul- 
ties arise with trying to construct a many-body variational wave function that 
^ I would take into account the locality of the flux-bubble interactions. In this 

C ' paper we look at a toy system, a mesonic molecule, in order to understand 

these difficulties. En route, a new variational wave function is proposed that 
v^ \ may have a sufficient impact on the old string-flip potential model results that 

;_( ■ the inclusion of perturbative effects may not be needed. 



1. Introduction 



For the past 30 years several attempts have been made, with little success, to 
describe nuclear matter in terms of its constituent quarks. The main difficulty 
is due to the non-perturbative nature of QCD. The only rigorous method for 
handling multi-quark systems to date is lattice QCD. However, this is very 
computationally intensive and given the magnitude of the problem it appears 
unlikely to be useful in the near future.^ As a result, more phenomenological 

^ Some very recent advancements have been made in the area of lattice QCD that have 
reduced computation time by several orders of magnitude. "Now what took hundreds of 
Cray Supercomputer hours can be done in only a few hours on a laptop computer" |1[. 
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means must be considered. 

A good phenomenological model should be able to at least reproduce qualita- 
tively all the overall bulk properties of nuclear matter. In particular, 

- nucleon gas at low densities with no van der Waals forces 

- nucleon binding at higher densities 

- nucleon swelling and saturation of nuclear forces with increasing density 

- quark gas at extremely high densities 

There are several models that attempt to reproduce these properties but 
none of them does so completely. In this paper, only the string-flip potential 
model p|-§l, will be considered. This model appears to be promising because 
it reproduces most of the aforementioned properties with the exception of 
nucleon binding: instead of producing a binding of about 8 MeV at nuclear 
density, one finds anti-binding of about 25 MeV. However, various efffects 
have been ignored, such as colour Coulomb and hyperfine forces, relativistic 
effects and many-quark clusters. In this paper we attempt to understand why, 
in terms of a toy mesonic molecular system P] . 



1.1. The String-Flip Potential Model 



The idea of string-flip potential models derives from lattice QCD and meson 
spectroscopy. A static potential derived from computations in lattice QCD 

V{r)r^ar-^^, (1) 

6 r 

is confirmed as a phenomenologically satisfactory potential model between 
quarks and antiquarks by fitting the experimental mesonic spectra. This is 
basically an interpolation between the long range non-perturbative {ar) and 
short range perturbative (— | ^) parts of the force between pairs of quarks 
(Fig. |l|). In the many-body case, the string- flip potential model ignores the 
short range part of the potential and considers an ensemble of quark- ant iquark 
pairs, qq, such that the total string length, J^^qq, is a minimum: i.e., 



V = CT Y. rqq. (2) 

inin{gg} 



This particular model has been used in an attempt to simulate nuclear matter. 
Although it has an obvious shortcoming, in that it is more applicable to a pion 
gas, it does surprisingly well at predicting some of the overall bulk properties 
of nuclear matter MM . 
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(a) (b) 

Fig. 1. The colour field lines between quarks collapse upon themselves, due to the self 
interacting nature of the gluons, to form a flux-tube-like structure. At long distances (a) 
the fields lines collapse to become almost string-like and at short distances (b) the fields 
lines expand to become almost QED-like. 

In a previous paper, we have extended this simple model to a more realistic one 
which involves triplets of quarks. Here the flux-tubes leaving each quark meets 
at a central vertex such that overall length, r, of "fiux-tubing" is minimized 
(c/. Fig. 0). The potential energy is simply ar 0. The effects of this are 
surprisingly small in going from SU{2) to SU{3): a nucleon gas still forms, 
and becomes a quark gas above a certain critical density. 

O 




Fig. 2. Flux-tube arrangements for the 3q cluster potentials. 



These models are motivated by results from lattice QCD, where variations 
are taken about minimal lattice field configurations between quarks. We will 
refer to these, "hnear potential models," as, ''SUe{N) models," where the 
SU{N) refers to the SUc{N) Yang-Mills gauge group and the subscript, "£," 
refers to, "linear potential," ^. In some models each quark has fixed colour 
to simplify the searches for minimal fiux-tube configurations as the system of 
quarks evolves. When necessary, for clarity, these models will be referred to 
as, ''SU'i,{N) models." 

This paper is divided into four sections: the first discusses the problems faced 
when attempts are made to extend the string-flip potential model to include 
perturbative QCD interactions for many-body systems. This is followed by 
a section on mesonic molecules, where the problems faced with the string- 
flip potential model and its extensions are examined in detail. In the next 
we examine a modified wave function, which appears to solve many of the 
difficulties. The final section reviews all of the findings and their possible 
consequences for modeling nuclear matter. 



2. The Flux-Bubble Model 



It appears from earlier work that the String-Flip potential model is incapable 
of producing nuclear binding without the addition of new interactions 0,|^ . 
However, it turns out to be surprisingly difficult to introduce these into a 
many-body model. A model was recently proposed, called the Flux-Bubble 
Model 0, which does allow the perturbative QCD effects to be included. 

The primary objective is to construct a model which combines both the non- 
perturbative (fiux-tubes) and perturbative (one-gluon exchange) aspects of 
QCD in a consistent fashion. In order to simplify this task only colour- Coulomb 
extensions to a linear potential model using SU{2) colour, i.e., SUi{2), will 
be considered. The conventional phenomenological potential. 



V{r) ^ ar- -— , 
6 r 



m 



mentioned in § |1.1| does not extend naturally to many-body systems, since the 
4/3 must be replaced by Ajj = —3/4, 1/4 for unlike and like colours respec- 
tively. This gives rise to a Van der Waals potential between colourless nucleons 
with a \ behaviour. The flux bubble modifies this to 



^(% - ro) 



if lij > ro 



V,j-~ 



^s'^ij 



if Tij < ro 




(3) 



Coulomb 



with as ~ 0.1 . The major difference is that the nonperturbative and pertur- 
bative regimes are completely separated in the latter, so that there are no 
Van Der Waals forces.|3 When the quarks are separated at a distance greater 
than ro the potential is purely linear and when they are inside this radius it 
is purely Coulomb. In effect, for distances less than ro, a "bubble" is formed 
in which the quarks are free to move around, in an asymptotically free fash- 
ion. In both distance regimes the net colour of the system is neutral, and 
phenomenologically the models are almost identical for a single qq system. 



This extension of the linear potential, although simple for a pair of quarks, 
becomes more complex when considering many pairs of quarks. In particular, 
it is difficult to construct a potential model when some set of quarks lie within 
ro, without forming a colour singlet so that they must be connected to more 
distant quarks by flux-tubes. An ansatz that satisfles these requirements is 
obtained by inserting virtual qq pairs across any of the intersection boundaries 

^ We have also considered a "smoothed" potential where the Coulomb term has an expo- 
nential cutoff of the form e ''o , and the linear term is turned on by a similar function. This 
makes a negligible difference to our results. 



formed by the flux-tubes with the bubbles. Now the segments of flux-tubes 
that lie outside the bubbles remain intact while the segments inside simply 
dissolve. Fig. ^ illustrates the dynamics of this model. 



Fig. 3. Consider configuration (a) of quarks, with r > ro, 
about to move to (b), s.t. two of them are within r < tq. 
Then the procedure is to draw a bubble of rg away from 
the two, (b), and to cut the flux-tubes at the boundary and 
insert virtual qq pairs, (c). Once the potential is computed 
the configuration is restored to (b) before the next move is 
made. 
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Cut here ~ ' Virtual qq pair insertion 



Note the insertion of the virtual qq pairs allows the construction of colourless 
objects. These are solely used as a tool to calculate the overall length of the 
flux-tube correctly, and not used in computing the Coulomb term however, as 
the fleld energy is already taken into account by the "real" quarks inside the 
bubbles. In general, once the bubbles have been determined, the flux-tubes 
must be reconflgured in order to minimize the linear part of the potential. 
Although the model is currently for SUe{2) it would be straightforward to 
extend it to a full SU£{3) model with all the one-gluon exchange phenomena. 

In the previous case for the SU^{3) model, given in ref. 0, we adopted a 
wave function with two "independent" parameters, p (density) and f3 (inverse 
correlation length); i.e., 

- E (P^nnT 

^a/3= e ^'^^^^^ n Mrp,)\, (4) 



Correlation <^ (3 Slater <^ p 

where a was flxed.[3 We were able to use a scaling trick to reduce the number 
of variation degrees of freedom to one, since p and /? varied parametrically 
with a single parameter, 9 , as 

(/3,pi/=^) — > C(^)(cos^,sin^) (5) 

where [C(^)] ~ fm~^. However, for the flux-bubble potential, this scaling 
•^ Here, a « 1.75, and a = 2 gives the harmonic oscillator case. 



transformation is broken. This extra degree of freedom greatly increases the 
computation time, since we must now perform a 2-dimensional minimization. 

We have, however, been able to show that the results of the computation 
are qualitatively similar to previous results : in other words, the flux bubble 
does not give rise to additional binding. However, to obtain satisfactory results 
requires a coarse 10 x 10 mesh of points, in p and /3, and hence 18hrs of CPU 
time on an 8 node farm. Clearly this procedure would be ridiculously slow if 
more parameters were to be added. Thus a rapid method of checking different 
wave functions, minimization procedures and various aspects of the string-flip 
and flux-bubble potential models is desirable. 



3. Mesonic Molecules 



A possible mini- laboratory is a mesonic- molecule [§|, Q2; consisting of two 
heavy quarks and two relatively light antiquarks: see Fig. ^.The quarks are 
assumed to be heavy so that the light antiquarks can move around freely 
without disturbing their positions: i.e., the Born-Oppenheimer approximation 
is valid. By varying the distance, R, between the heavy quarks a mesonic- 
molecular potential, U{R) , can be computed. 




q 

Fig. 4. A mesonic-molecule, Q2 , with two heavy quarks and two light antiquarks. 
The Schrodinger equation that describes the effective potential, U{R), is given 

byi 

^ Ev; + V^|^ = t/(i?)^; (6) 



2"^g q 



we shall assume m^ ~ 330 MeV. The potential V describes the many-body 
nature of the four quark system, and is therefore model dependent. This equa- 
tion can be solved variationally for U{R), at fixed values of R, by guessing at 



the form of the wave function, \l/, and minimizing 

U = f + V (7) 

with respect to the parameters in \1/. 

The expectation values T and V are found by using the Metropohs algorithm 
PJTO|], and the optimal parameters for \Ef are found by using the distributed 
minimization algorithm [p|,pA| on U{R). This was developed to handle the 
problem of reducing CPU overhead. It involves the usage of several CPU^s to 
perform interlaced parabolic searches for the minima of a given surface. For 
the Q2 system this algorithm yields 

^CPU - 0(£) , (8) 

for p parameters and m computers: similar results are expected for nuclear 
matter calculations. 

The Q2 system provides a good way of checking potential models for the 
possibility of nuclear (mesonic) binding. Because of its simplicity, it allows for 
investigating various wave functions and minimization schemes with minimal 
computational effort. 

A final useful feature is that the integration can be reduced to a 4-D one if 
Coulomb interactions are ignored, and a 5-D one if they are included. This 
allows us to avoid the Metropolis procedure entirely, and use multi-dimensional 
Gaussian methods, which we use a check on the Monte-Carlo methods 

3.1. Sensitivity to the Variational Parameters 



When the SUe{2) string- flip potential model was investigated in ref. 0] the 
parameter a in the variational wave function, given by Eq. ^, was fixed by 
requiring that it minimize the total energy at zero density. Since this could 
be done analytically it allowed for a reduction in the number of variational 
parameters used in the Monte Carlo. It was assumed that constraining a would 
have very little effect on the physics as a function of density since the results 
only vary by about 1% for 1.5 < a < 2.1 at zero density. The validity of this 
claim can now be checked more thoroughly by using the Q2 mini-laboratory. 

Fig. H shows a plot of the Monte Carlo results, obtained via Eqs. ^ P, and |l^, 
for f7(-R) where a is allowed to vary, a = 2.00 , and a = 1.74 . The a = 2.00 
value was used in the old SUi{2) model |^, and a ~ 1.74 was the value that 
minimized U{R) at infinite separation. The values of U{R) at the end points 
of the curves, from R = fm and out to i? = 5 fm, were checked against the 
analytic solution given by Eq. |A.8 . 
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Fig. 5. t/(i?) where a is, allowed to vary, fixed at 2, and fixed at 1.74. 

We find that when a is allowed to vary, a and /? become highly correlated. 
Regardless, looking at the maximum fluctuations about the central values of 
the parameters a, /3, and d{R) we have § a ^ (1.88 ± 0.13), /5 ^ (1.32 ± 
0.05) /m, and t/(i?) ■^ (1710 ±4) MeV , which is certainly less then a 1% effect 
in the total energy [|]. However the relative energy U{K) — U{oo) depends 
critically on the variational procedure. 

Surprisingly U{R)\a=2 gives a much deeper well than if a was left as a free 
parameter: i.e., D ^ C^(8) MeV [^]. For certain fixed values of a , for example 
a ~ 1.74 in Fig. |^, the potential gives a slight short range repulsion, which 
becomes quite dramatic for a ~ 1.0. This is, of course, totally unphysical, 
since U{R) — f/(oo) < 0. We conclude from this that it is unsafe to fix a , 
since the well depth can vary by nearly a factor of 3. However, it is fairly safe 
to assume that the results of past papers 0J^ will not change significantly 
if a is allowed to vary. For the rest of the paper, we consider simultaneous 
minimizations of a and /3 . 



3.2. A General Survey of Extensions to SUi(2) 

In this section the effects of extending the old SUe{2) 0] model to include fiux- 
bubbles, with and without fixed colour {i.e., SUe{2) and SU'^{2) respectively), 
will be investigated in the context of the mesonic-molecular system, Q2 ■ For 
the old SUe{2) model the variational wave function was assumed to be of the 
form |[, 

^„^ = e ™>^{Qq} ^ (9) 

where a and f3 are variational parameters, and rqg is the distance between a 
given Qq pair. The summation in the exponent of the wave function is over 
the set of Qq pairs that requires the least amount of flux-tubing: i.e., 

^ = ^ E rQq = (rmm{ ^ rgg-} , (10) 

min{Qq} {Qq} 

or more explicitly 

V = amin{rQ,g-, + rg^ga, rQ^q2 + rg^qA , (H) 

cf. Eq. (0). Therefore, the kinetic energy is simply 

T=£-{ T. Hl-iPrQ,r) + l]r-Qf). (12) 

^'"9 ram{Qq} 



Now, if we assume for the moment that the colour is affixed to the quarks 
{i.e., SU[{2)), then the most general fiux-bubble potential is of the form 

V = a Y. ('^Qq- - '^o) 0{rQq - tq) + a, ^ Kpj ^(^0 - r^^p, ) , 

min{Qg} i<j VPiP] '^O / 

(13) 
with particle index p^ e {Qi, qj\i,j = 1, 2} , such that k = 1,2,3,4 , and SUc{2) 
colour factor 

A,«4-: ''""''"*'*' . (14) 

[ 4 if PiPj £ {qq, QQ} 

This potential can be rewritten into the more enlightening form 



mm{qiqj} 



^ (rq.q, - '^o) 0{rg^g^ _ ro) - - a. e{ro - r^^g/^ 



min{qiqj}. ^ ^ ■< 



qiqjemin{qiqj}_^_ 



+I«^ ^_ f r^ - 7) ^(-0 - r,,J , (15) 

— ^~r — ^T \' Qiqi '0/ 



where g^e {Q, q} ■ The sets minjgjgj} and min{gjgj}_ contain the attractive 
and repulsive quark pairings, respectively, in the complement (min{gjgj}) of 
the set min{gjgj} : i.e., 



mm{qiqj} = mm{qiqj}^ U mm{qiqj}_ .Q 

The two terms inside the square brackets of Eq. |T5| represents a linear potential 
followed by its colour- Coulomb extension: cf. Eq. (0). The remaining terms 
represent the rest of the colour- Coulomb interactions, which contain both 
attractive and repulsive bits. 

Fig. ^ shows the results of the Monte Carlo for the linear potential and its 
various variants, including the flux-bubble potential. 

We have also considered relativistic effects, via the relatvistic Dirac-Hartree- 
Fock method(e.5i., [0), assuming the fiux-tube interaction is a Lorentz scalar. 

E [^l + im, + arQ,f]^ = U'iR)^. (16) 

minjQij} 

This slightly exacerbates the problem, in that the potential becomes even more 
more shallow. We note in passing that it is not possible to solve the physically 
interesting case of massless quarks interacting via a Coulomb potential in 
addition to the linear term, since it introduces an effective 1/r^ potential in 
Eq. 



To compare these potentials with more conventional methods, we would like 
to investigate whether they can bind heavy quarks into molecular systems, 
and to what extent they resemble conventional inter-nucleon potentials. The 
details of this analysis can be found in appendix 0. A reduced heavy quark 
mass of 

I^Q> (53.7 ± 1.9) GeV (17) 

is required to obtain binding for these potentials. This should not be too 
surprising as the potential wells are very shallow, 

i) ^ (2.986 ± 0.030) MeV^ . (18) 

To find whether the string-flip potential models actually mimic pion exchange, 
the asymptotic parts of the U{R) 's are fitted to a Yukawa potential. The results 



.g., if inm{qiqj} = {Qiq2, Q2<7i} then mm{q.iqj} = {Qiqi, Q2q2, Q1Q2, 9192}, which 



imphes mm{qiqj}_^ = {QiQi, Q2q2} and m.[n{qiqj}_ = {Q1Q2, 9192}- 
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Fig. 6. The binding energy, Eb{R) (= U{R) — Uoo s.t. C/qo = U{oo)), for the Q2 system as 
function of heavy quark separation for a linear (Eq. (Eof)), a linear-plus-Coulomb (square 
brackets of Eq. (p^)), and a flux-bubble (Eq. (|l5|)) potential model. Also included, for 
contrast, are the semi-relativistic results discussed in § ^. The curves in these plots have been 
parameterized by u(r) - Mo{e"^''''(''~'^«''' -2Ke"''(''"'^")°}(r--e)'' (c/. Morse potential §,|l|). 

of the analysis yielded an effective exchange mass of 

mex ~ (636 ± 41) Me\/ , (19) 

which is about 4.5 times too big, so the potential is much too short-range. 

Finally, the flux-bubble potential is extended to allow the colour to move 
around. The potential is similar to that of Eq. |1^ except now the particle 
indices, p^, carry colour degrees of freedom: 



qu e {Ri,Bi,ri,bi\t = 1,2} C ^^7^(2) 
qkS {R,,B,,f,,h\t = 1,2} C SUc{2) 



11 



I.e., 

fe~f,...,etc. ^ SU,{2) ~ SUc{2), 

where the capital case letters represent the heavy quarks, the lower case letters 
represent the light quarks, the letters r and R represent the red quarks, and 
the letters h and B represent the blue quarks. Therefore, the minjQg} in 
Eq. |l^ is determined by Eq. [10| such that colour is no longer fixed to a given 
quark. Fig. |^ shows the results of the Monte Carlo for this potential. 
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Fig. 7. The Q2 binding energy curve using the fuU flux-bubble interaction, i.e., in which 
the colour is allowed to move around. This curve was parameterized using u{r) given in 
Fig. |. 



Obviously the result of this is totally different from the fixed colour model. 
For R <1 fm, the potential just reproduces the potential between the heavy 
quarks: a linear fit to this region (see the insert in Fig. |^ yields a slope of 
(909.8 ± l.Q)MeV/ fm\ Near the origin the potential behaves as — -, which 
is due to the Coulomb attraction between the two heavy quarks. Beyond 
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R = 1 fm there is a barrier and beyond this no more structure. This apparently 
bizarre result can be understood as follows: if two mesons, each containing a 
heavy quark and a light anti-quark, are brought together from infinity they ini- 
tially feel a repulsive force. However, it then becomes energetically preferable 
for the two mesons to dissociate into one meson containing two heavy quarks 
and another containing two light anti-quarks. This arises since the model has 
a Pauli-Gursey symmetry |]n[: quarks and antiquarks have the same repre- 



sentations. Although potentially realistic as a description of meson-anti-meson 
interaction, the SUc{2) with moving colour does not make a satisfactory model 
for nuclear matter. 

From this discussion, it is clear that the extensions to the string-flip potential 
model to include the colour- Coulomb interactions have essentially no effect on 
the Q2 potential, and certainly will not give rise to nuclear binding. Further- 
more, the string-flip potential does not give a long-range interaction similar 
to pion exchange. When the colour was allowed to move around a highly un- 
physical situation occurs which suggested that there was perhaps a problem 
with using SUc{2) or with the variational wave function itself — perhaps even 
both. In the next section we investigate an alternative wave-function. 



4. A New Wave Function 



The two previous sections have shown that modifying the potential or mod- 
ifying the variational procedure produces a relatively small effect on the Q2 
interaction. It therefore seems plausible that the fundamental problem lies in 
our choice of the wave-function. A hint is to note the similarity between the 
Q2 mesonic-molecular system and H2 molecular system. The key reason for 
the molecular binding is the screening effect caused by the electrons which 
are for the most part localized in between the protons. This localization is 
achieved by using a variational wave function that is a superposition of the 
direct product of two ground state hydrogen atoms, [|,|l^,|l^ 



^ ~ e + e . (20) 

The effect is that of a "bond-centred" wave function. Although the Q2 system 
is far removed from its H2 cousin from a dynamical point of view and the mo- 
tivations for achieving localization are quite different, it would seem plausible 
to use a similar ansatz: 

^a,f3 = e +e . (21) 

If qiQi and q2Q2 represent two separate mesons then the first term represents 
the internal meson interactions while the second term represents the exter- 



13 



nal meson interactions. Notice that the external interactions shut off as the 
separation, R , between the two heavy quarks becomes large, 



lim '^{R) = e 

R—too 



-nr^^Q.+^lQ.) 



f22) 



which is the desired property. Furthermore, when the light quarks are close in 
space, ^(i?) is considerably enhanced. 



Using Eq. |^ and Eq. ^ the kinetic energy contribution is more complex for 
this wave-function, and now becomes ||,|,|l^] 



where 



and 



T = 2T_s- F\ 



(23) 



T-s = --U^lln^) 
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(24) 
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-/^"('"flQl + '^IQi + '"51Q2 + '^Iq2)\ /^p.^ 

xe ) . (25) 

It serves as a good check of the computation that in the large R hmit the 
kinetic energy reduces to the kinetic energy for the old wave function in the 
same limit: 



"1 
I 2 



2^^\EKi-(/5-.-.Q.r) + i]CQ!];^ (26) 

c/. Eq. 0. Also direct evaluation of the RHS leads to (assuming ^ is properly 
normalized) 

limf(/^) = ^^/5^ (27) 

which is just the kinetic term for the analytic solution given by Eq. |A.1| . 

The Monte Carlo computations that were done, using the pseudo-hydrogen 



wave function, '^ h [i.e., Eq. ^ in § |3.2| , have been repeated here for the 



pseudo-hydrogen-molecular wave function, '^h2 (^-c, Eq. |2l]), and are shown 
in Figs. H and |. 

It can be immediately seen, that there is a dramatic contrast between the 
figures for ^h and ^h2 ■ The depth of the pseudo-Morse potential has in- 
creased from around 5 MeV to around 80 MeV. This is sufficient to bind 
quark molecules. From the plots in Fig. ^ we find that heavy quarks with 

^iQ > (660 ± 24) MeV . (28) 

such as c and h would form bound quark molecules. We can fit a Yukawa 
potential to the tail of the potential to give an effective Yukawa mass 

mex ~ (575 ± 32) MeV , (29) 

via table [A.3| , which is still high compared to the pion mass. Also included in 
Fig||, are plots for a fixed at 2 and 1.74. The subtle effects seen with the old 
wave function are simply overwhelmed by the depth of the potential: in fact, 
all of the curves are independent of the details of the potential. Therefore, 
any new model of nuclear matter that incorporates '^ H2 should run into no 
difficulties by fixing a. 

The final figure. Fig. p, of the fiux-bubble model with moving colour is quite 
intriguing. The anomalies in Fig. ^ have disappeared; the light quarks have 



15 



-10 



-20 



-30 



-40 



-50 - 



-60 - 



-70 - 



-80 - 




-90 



*-- 


Linear 


--0-- 


Linear + Coulomb 


^&- 


- Flux-Bubble 


- - -A- 


- a = 2.00 


#- 


a = 1.74 


0- 


Composite 



0.25 0.5 0.75 1 



1.25 



1.5 



1.75 



.25 2.5 

R (fm) 



Fig. 8. The binding energy, Eb{R), for the Q2 system as function of heavy quark separation 
for a linear, a linear-plus-Coulomb, and a flux-bubble potential model. Also shown, are the 
curves for a = 2.00 and 1.74 , and the composite model discussed in § g[ The curves in these 
plots have been parameterized by u{r) given in Fig g. 

not drifted away as an isolated pair to leave a linear potential between the 
heavy quarks. Naively, this seems to suggest that the problem was with the 
wave function and not SUc{2). However, this is not quite the case, since the old 
wave function gives an interior well depth which is twice as deep. Therefore, it 
is still energetically more preferable for the Q2 system to dissociate into two 
isolated mesons; one with two light quarks and the other with two heavy ones. 
The interior of the well can be fitted to a Coulomb potential of the form. 



Vc{r) = - 




(30) 



(see insert in Fig. D, with a ^ (135.5 ± 2.9) MeV fm, (3 ^ (2.844 ± 0.048) 
fm~^, Voo ~ (68.3 ±2.6) MeV. The term in the brackets is included to mimic 
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Fig. 9. The Q2 binding energy curve using the full flux-bubble interaction, i.e., in which 
the colour is allowed to move around. This curve was parameterized using u{r) given in 
Fig. |. 



the overlap between the charge distributions of two mesonic systems. In terms 
of as, the hyperfine constant for this region of the potential is 



ag^ ^ (6.86±0.15)as. 



(31) 



The interior part of the potential is quite deep and bottoms out at C(— 270) 
MeV, at which point the —ag/r term for the heavy quarks kicks in {i.e., for 
R ^ 0.1 fm). The exterior part of the potential fits to a Yukawa potential 
with m,, ^ (612 ± 32) MeV (wa table |A^ ). 
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5. Discussion 



Various aspects of model building for nuclear matter have been examined in 
the context of the Q2 system. The ramifications of these investigations will 
now be discussed. The most important result is the new variational wave func- 
tion, Eq. |2l|. This new wave function makes a massive change in the depth of 
the Q2 potential well, increasing it by a factor of 27, which is deep enough to 
bind relatively light quarks: i.e., mq>jO{ms). The wave function also fulfills 
the requirement of handling local flux-bubble interactions, which becomes ap- 
parent when looking at the "before" and "after" pictures (of moving colour) 
in Figs. and P, respectively. Therefore, in SUi{2) for a many quark system 
this would suggest the following ansatz: 



vl>~Perm|^H(rp.p,)| U |<l>(r,J| (32) 

colour 



colour 

where 

^Hir,^,,) = e ' ""'" , (33) 

Perm|\l'//(rp^p^-)| is a totally symmetric pseudo-hydrogen wave function and 
|$(rp^,)| is a totally antisymmetric Slater wave function. 

For the full three-quark system a similar wave function would apply. This does 
not necessarily mean that this wave function would lead to nuclear binding, 
but, given the order of magnitude of increase in the Q2 well depth it would 
seem quite plausible that it might be a strong enough effect to produce a 
shallow well in the nuclear-binding-energy curve. A simple test would be to 
consider SU(,{2) with just a string-flip potential, with a flxed, in which case 
scaling is restored and the Monte Carlo becomes quite straightforward to do. 

When the SUi{2) flux-bubble model with moving colour was considered, the 
results showed that the Q2 system dissociated into one light and one heavy 
meson. However, this model is not physical: rather we should consider the 
heavy quarks as a composite of two light quarks and use 5'f/c(3) instead, so that 
a flux-tube cannot form between the two heavy quarks. For this "composite" 
model, the Ap.p^. 's of potential Eq. |13| become 

\ if PjPj e {66, 55} 

I if ViVj ^ {^G, Bg} 
-liip,p,e{bB,gG} , (34) 

-| a piPje{bg,BG} 
-| a PiPjE min{g,.gj 

such that rg r^ B and rb ^^ G for the composite states, bb ~ gg and BB ~ GG. 



^PiPj 



If the heavy quarks are considered to be a composite of two hght u quarks 
then the effective mass would be about 600 MeV, which is a httle too small to 
produce binding in this model, but for masses this low the adiabatic approx- 
imation would no longer be valid. Therefore, an interesting possibility would 
be to consider a many-body SUe{2) or SUe{3) flux-bubble model in which 
there is an imbalance between the quark and anti-quark masses. An inter- 
esting side effect of the improved wave-function is that it will automatically 
produce nucleon deformations, which, it has been argued, produce the contact 
interactions responsible for nuclear binding W^. 



Further enhancements are expected by going to the full qqq nucleon model, fol- 
lowed by including flux-bubble interactions. In particular, the colour-hj^erfine 
interaction which is shown by linked cluster expansion models to play a sig- 



nificant role in nuclear processes [|I3- The fiux-bubble model proved quite 



successful at combining colour-Coulomb interactions with fiux-tube interac- 
tions. Although these interactions had very little effect on the Q2 system, it 
was useful in demonstrating that extending the flux-tube model to include 
local perturbative QCD interactions can be done. This is important, since the 
nuclear hard-core potential is usually believed to have its origin in hyperfine 
interactions [3,^|T7| . Therefore, it would prove most interesting to investigate 
the effect of adding more perturbative interactions to the Q2 system. It also 
appears that relativistic effects are unimportant. With the addition of SUc{3) 
this would lead to a more realistic model of mesonic molecules which perhaps 
could be tested in the laboratory. 

The Q2 system has proven to be a very useful aid for trying to sort out the 
complexities of model building for nuclear matter. The details of the mechan- 
ics, from wave functions to dynamics to practical computing methods, of the 
flux-bubble model have now been thoroughly investigated. It appears that 
the flux-bubble model may prove to be very successful, not only for modeling 
nuclear matter but also for modeling mesonic molecules as well. 
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Appendix A. Analysis 



This appendix contains a summary of the analysis done for the models given 
in sections ^^ through |^ of this paper. 



All of the results of these sections make use of the general potential Eq. |T^. 
For lf{R), as defined by Eq. |^, some analytical results were obtainable where 
the Q2 system effectively dissociates into two isolated meson systems. For all 
of the plots this occurs at U{oo), and for the linear and linear-plus-Coulomb 
plots this also occurs at f/(0). In these regions, the analytic expression is given 

by (c/. i) 



where 



such that 



E 



free 



2/i 






^o(a,/5)ro 



--as \gc{a,^)^ 



(7ja, /3) = [1 - P(4/a, 2(/5ror ) ] ^ J«) 
^o(a,/3) = l-P(3/a,2(/5ro)"), 






^c(a) 



^-2 
a^2" r(2 + l/a) 

r(3/a) 

r(4/«) 

2°r(3/a) 

r(2/a) 
_i ' 

2 "r(3/a) 



(A.l) 



(A.2) 
(A.3) 
(A.4) 



(A.5) 
(A.6) 

(A.7) 



P(a, 2;) = 1 — r(a, z)/r(a), and T{a,z) is the incomplete gamma function 
pl| , |22[] . In the limits as tq — > and tq -^ 00 Eq. ( |A.1| ) reduces to solutions for 
the purely linear, 



Er * 

J-ifree ^ ^ 

and purely Coulomb, 



E 



free 






' '-^ ^' -\asgc{a) ^y 



(A.8) 



(A.9) 
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cases, respectively. Table [A.l| gives a summary of the Monte Carlo vs. an- 
alytic results. The results for the linear and linear-plus-Coulomb potentials, 

Table A.l 
Monte Carlo (MC) vs. analytic results for Figs. |[ |[ |and|. 



Method 




Fig. 


R{fm) 


Efree {MeV) 


a 


(3 (/m-i) 




Analytic 


Eq. A. 8 


n.a. 


O/oo 


1709.61 


1.75 


1.37 



MC 



Linear 



n 





1709.424 ± 0.038 


1.74 


1.37 




5 


1709.492 ± 0.044 


1.72 


1.38 


1 





1709.647 ± 0.032 


1.77 


1.36 




5 


1709.485 ± 0.032 


1.75 


1.36 



Analytic Eq. A.l 



n.a. 



O/oo 



1527.07 



1.74 



1.37 



MC Linear 

+ 
Coulomb 



1527.171 ±0.029 1.78 1.36 

5 1527.171 ±0.029 1.78 1.36 

1527.884 ± 0.041 1.74 1.34 

5 1527.884 ± 0.041 1.74 1.34 



MC 



Flux-Bubble 



1 


5 


1526.974 ± 0.032 


1.74 


1.36 





5 


1527.884 ± 0.041 


1.74 


1.34 


1 


5 


1527.884 ± 0.041 


1.74 


1.34 



1527.173 ±0.062 1.68 



1.40 



at i? = fm and R = 5 fm {i.e., ~ oo), were checked against Eqs. K^ 
and |A.1| , respectively. The fiux-bubble case, at i? = 5 fm, was verified using 
Eq. [A.l| . The minima of the analytic expressions, which were used to verify 
the aforementioned models, were found via the FindMinimum[. . .] routine in 



Mathematica 21 



The binding energy of the Q2 system can be estimated by doing a local 



parabolic fit ||T2[ about the minimum of U{R): i.e., by fitting 

y{r) = C{r - rof - V (A.IO) 

such that C = ^Quj'^/[2{hc)'^], where fiq is the reduced mass of the heavy 
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quarks. Therefore, the binding energy is simply 

E^, = iu+^)hco-V, (A.ll) 

which imphes 

/iQ > /imin = -^^ (A. 12) 



in order to obtain binding. Table |A.2| shows the results for the parabolic fits 
to the U{Rys in Fig. |. 

Table A. 2 

Results for local parabolic fitting about the minima of the plots in Figs. and 0. 

Model Fig. C{MeV/fm^) ro{fm) V{MeV) /imin (GeF) 

Linear | 29.8 ±1.5 0.4174 ± 0.0056 3.268 ± 0.058 54.3 ± 3.3 

I 234.8 ±8.0 0.5288 ±0.0057 82.55 ± 0.49 671 ± 24 

Linear ± | 28.7 ±1.4 0.4134 ± 0.0045 3.098 ± 0.060 58.2 ± 3.6 

Coulomb I 230.3 ±8.0 0.5833 ± 0.0058 82.60 ± 0.49 657 ± 24 

Flux-Bubble | 20.0 ± 1.0 0.4614 ± 0.0034 2.783 ± 0.042 50.3 ± 2.9 

226.8 ±8.0 0.5840 ± 0.0059 82.24 ± 0.49 653 ± 24 



The asymptotic parts of the U{R)^s were fitted to a Yukawa potential of the 
form 

VY{r) = -2f . (A.13) 

r 

Table |A.3| summarizes the Yukawa fits for all of the models. 
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Table A.3 

Results for asymptotic Yukawa fits of the plots in Figs. 0, 0, and S, where jtIox = hen. 

Model Fig. f {MeV fm) fJ.{fm-'^) m^^^MeV) 

Linear | 2.49 ± 0.20 3.23 ± 0.21 637 ± 41 

I 27.3 ±2.7 2.97 ±0.16 586 ± 32 

Linear ± Coulomb | 2.0 ±2.2 3.0 ±3.5 590 ± 690 

I 25.9 ±2.7 2.89 ±0.16 570 ± 32 

Flux-Bubble | 1.8 ±2.3 2.4 ±3.5 470 ± 690 

I 25.9 ±2.7 2.89 ±0.16 570 ± 32 

I 30.3 ±2.1 3.10 ±0.11 612 ±32 
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